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1  Introduction 


The  pursuit  of  enhanced  aircraft  performance  is  an  ongoing  venture  for  the  Air 
Force,  NASA,  and  the  aviation  in  lustry.  Their  efforts,  and  the  efforts  of  others, 
are  yielding  numerous  new  technologies  which  are  rejuvenating  many  previously 
abandoned  progrms.  One  such  program  is  the  design  and  development  of  a  profitable 
High  Speed  Transport  (HST).  future  HST’s.  It  ha'  been  concluded  that  for  such  on 
aircraft  to  be  economically  viable,  significant  performance  improvements  over  those 
offered  by  current  technologies  must  be  made. 

One  of  the  most  promising  prospects  for  performance  gains  on  HST’s  is  the 
skin  friction  drag  reduction  associated  with  obtaining  laminar  flow  over  the  wing.  It 
has  been  predicted  that  laminar  flow  control  (LFC)  would  provide  more  dramatic 
improvements  than  any  of  the  other  possible  technolog’  cal  advances  in  aerodynamics 
(reference  1).  The  most  encouraging  means  of  doing  this  is  through  the  use  of 
Hybrid  Laminar  Flow  Control  (HLFC)  systems.  In  such  a  system,  the  wing  is 
designed  to  produce  a  pressure  distribution  conducive  to  natural  laminar  flow  (NLF) 
over  its  midsection,  while  in  the  leading  edge  area  the  flow  is  stabilized  by  suction 
through  perforated  panels  or  narrow  slots.  HLFC  research  has  already  established 
an  advantage  in  drag  reauction  for  subsonic  aircra*.,  and  if  significant  amounts  of 
laminar  flow  can  be  maintained  on  HST’s,  where  skin  friction  typically  accounts  for 
more  than  oue-third  of  the  cruise  drag,  it  could  greatly  enhance  range  and  payload 
capabilities.  The  associated  reduction  in  aerodynamic  heating  would  also  significantly 
reduce  the  structural  weight,  minimize  the  required  thermal  managemenr,  and  prolong 
the  material  life.  For  these  reasons,  supersonic  laminar  flow  control  has  been  targeted 
by  NASA  as  a  crucial  research  area  in  the  development  of  future  High  Speed  Civil 
Transport’s  (reference  2). 

There  are,  however,  several  adverse  effects  of  designing  aircraft  to  accommo¬ 
date  hybrid  laminar  flow  control  systems.  These  include  the  added  weight  and  power 
requirements  of  the  air-suction  system,  the  possible  maintenance  difficulties  associated 
with  the  porous  suction  walls,  and  the  probable  aerodynamic  performance  penalty 
of  sustaining  an  upper  surface  pressure  distribution  which  is  amenable  to  HXFC.  In 
this  paper,  the  latter  is  addressed.  A  wing  which  was  cambered  solely  to  possess  a 
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pressure  distribution  favorable  to  IILFC  would  theoretically  have  more  drag-due- to- 
lift  than  would  be  inherently  produced  by  a  wing  that  was  cambered  for  minimum 
drag-due-to-lift.  disregarding  any  boundary  layer  considerations.  The  intent  of  this 
study  is  to  quantify  tliis  difference  and  determine  its  likely  impact  on  the  true  benefits 
of  HLFC  wing  designs. 

It  was  determined  that  a  parametric  study  of  supersonic  arrow  wings  would 
be  useful  in  assessing  these  effects.  The  influence  of  tar  cr  ratio  and  aspect  ratio 
were  examined  by  including  several  planfcrms  of  varying  shape  in  the  study.  The 
investigation  was  conducted  at  a  Mach  number  of  2.0.  since  the  effect  of  camber  on 
minimum  drag-due-to-lift  decreases  at  higher  Mach  numbers.  Thickness  profiles  of 
the  wings  are  based  on  NACA  64-series  airfoils.  Two  different  design  philosophies, 
a  loading  optimization  method  and  a  method  to  design  HLFC  'wings,  were  used 
fo  camber  each  planform.  The  first,  more  conventional,  design  method  involves 
optimizing  wing  camber  and  twist  for  minimum  drag-due-to-lift  at  the  design  lift 
coefficient  with  no  consideration  for  boundary  layer  type.  This  method  has  been 
used  often  in  the  preliminary  design  process,  recent  examples  of  which  are  described 
in  references  3-5.  The  second  design  method  involved  cambering  the  wing  to  produce 
a  given  upper  surface  pressure  distribution  that  had  been  determined  to  be  desirable 
for  HLFC  wings.  To  achieve  this,  the  same  linear  theory  design  program  which 
was  used  to  design  the  wings  for  minimum  drag-due-to-lift  was  implemented,  since 
programs  which  use  higher  order  methods  to  design  a  three  dimensional  wing  for  a 
given  pressure  distribution  are  still  in  the  development  stages.  This  program  has  been 
found  to  be  adequate  in  many  previous  preliminary  designs,  while  providing  the  short 
computational  turn-around  times  required  to  conduct  a  parametric  study  on  several 
planforms  in  a  limited  time  period.  The  current  research  has  resulted  in  a  process  by 
which  the  program’s  optimization  routine  can  be  used  to  camber  a  wing  such  that 
it  produces  a  spanwise  lift  distribution  which  will  help  to  minimize  drag-due-to-lift, 
while  maintaining  any  given  upper  surface  pressure  profile  in  the  streamwisc  direction, 
such  as  that  desired  for  HLFC  wings. 
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2  Prelimenary  Evaluation  of 
Sr.personic  Wings 


For  the  evaluation  of  aerodynamic  performance,  the  force  estimates  will  be  built  up 
through  superposition  as  shown  in  Figure  2.1.  For  this  study,  we  will  consider  the 
drag  components  of  skin  friction  drag,  zero-lifi  wave  drag,  and  drag-due- to- lift,  on  the 
wing  only.  The  drag-due-to-lift  will  include  finite  wing  and  leading  edge  vortex  thrust 
effects.  All  of  the  methods  used  have  been  correlated  with  experimental  data  in  several 
NASA  studies  including  that  of  reference  6,  which  compares  three  computational 
methods  to  wind  tunnel  test  data. 


2.1  Zero-Lift  Wave  Drag 

The  wave  drag  of  a  configuration  is  typically  calculated  in  one  of  two  ways;  either 
by  the  far-field  (supersonic  area  rule)  method  or  by  the  near-field  (surface  pressure 
integration)  method.  Both  wave  drag  estimation  methods  have  been  used  to  evaluate 
the  wings  presented,  and  the  results  are  compared  for  equivalence. 

2.1.1  Far- Field  Wave  Drag 

The  zero  lift  far-field  wave  drag  is  computed  using  the  supersonic  area  rule  outlined  in 
reference  7.  The  supersonic  area  rule  is  a  derivative  of  the  transonic  area  rule,  which 
asserts  that  the  wave  drag  of  an  aircraft  is  essentially  the  same  as  the  wave  drag  of 
an  equivalent  body  of  revolution  having  the  same  cross-sectional  area  distribution. 
However,  in  the  supersonic  area  rule,  several  equivalent  bodies  of  revolution  are 
produced  by  passing  a  series  of  para’lel  cutting  planes  through  the  configuration, 
inclined  at  the  Mach  angle  a  to  the  aircraft  axis  (Figure  2.2).  Each  equivalent  body 
is  calculated  at  a  constant  azimuth  (roll)  angle  9.  Then  the  area  of  the  equivalent- 
body  at  each  longitudinal  station  is  the  projection  of  the  area  intercepted  by  the 
Mach  plane  onto  a  plane  normal  to  the  aircraft  axis.  The  wave  drag  for  each  of  these 
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equivalent  bodies  at  given  free  stream  conditions  is  then  calculated  using  the  Von 
Karman  slcnder-body  formula. 


r:1 

0(0)  =  ZLZisa.  j  J  4"(Xl)  a»(x2)  log  I  .r,  -  *2  j  dxx  dx7  (2.1) 

0  0 

where  i  is  the  length  of  the  aircraft,  A"  is  the  second  derivative  of  the  body  area 
distribution,  and  xx  and  x7  are  the  axial  locations  of  two  adjacent  Mach  planes  which 
enclose  the  volume  clement.  The  total  wave  drag  of  the  entire  configuration  is  then 
computed  as  the  integrated  average  of  all  of  the  equivalent-body  wave  drags. 


o 


In  the  program  used  to  perform  the  calculations,  equations  2.1  and  2.2  are  integrated 
numerically  over  a  specified  number  of  intervals  along  the  ar-axis  and  in  roll  angle  0. 
The  program  also  assumes  the  aircraft  to  be  laterally  summetric,  thus  calculations 
only  need  to  be  made  for  0  <  0  <  r.  For  this  study,  60  intervals  in  the  x  direction 
and  16  intervals  in  6  (over  a  range  of  it)  are  used. 


2.1.2  Near-Field  Wave  Drag 


The  near-field  method  of  computing  wave  drag  is  based  on  the  integration  of  zero-lift 
thickness  pressures  over  the  configuration  surface.  The  velocity  potential  over  the 
wing  is  calculated  using  the  familiar  linear  theory  equation 


1  f  f _ A  dg  d( _ 

*  •'  l  \/(x  -  <)’  -  Wy  -  v)2 


where 

6(x,y)  =  velocity  potential  at  wing  field  point  (x.y) 

A  =  surface  slope  at  integration  point 
3  =  vA/2  —  1  (supersonic) 

f  —  x  variable  of  integration 

t)  -  y  variable  of  integration 

r  =  interval  of  integration  defined  by  the  surface  of 
the  .  ing  within  the  Mach  cone  from  (x,y). 


(2.3) 
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The  wing  tliickness  pressure  coefficient  is  calculated  using  the  relationship 

P  -  Pi c 


Cp{x.y)  - 


(ho 


_ 

tfoo  dx 


(2.4) 


where 

/;  =  local  pressure  at  (x.y) 

Poo  =  free  stream  static  pressure 
u  =  x  perturbation  velocity 
£/*,  =  free  stream  velocity 
ijoc  =  free  stream  dynamic  pressure 

The  computer  programs  that  were  used  for  the  supersonic  analysis  have  this 
velocity  potential  integration  approximated  by  a  numerical  summation  on  a  system 
of  rectangular  elements  known  as  a  "Mach  box  '  grid  system  (fig  2.3).  A  field  point 
is  defined  by  its  coordinates  {x.$y)  and  an  influencing  element  is  defined  by  the 
coordinates  (£,,dp).  Note  that  only  the  elements  ahead  of  an  inverse  Mach  line 
originating  at  the  center-rear  of  the  field  point  (the  shaded  area)  are  included  in 
the  summation.  In  the  calculation  of  surface  pressures  due  to  thickness,  the  wing  is 
completely  defined  by  its  thickness  distribution  (to  define  the  surface  slopes  A).  The 
near-field  wave  drag  is  then  calculated  by  integrating  the  pressure  components  in  the 
axial  direction. 


2.2  Drag-Due-to-Lift 

For  a  first  order  analysis,  the  drag-due- to- lift  is  also  computed  using  linear  aerody¬ 
namic  theory.  The  wing  is  evaluated  using  linear  theory  solutions  (equations  2.3  and 
2.4)  on  the  same  type  of  “Mach  box”  grid  used  in  the  near-field  wave  drag  solution. 
However,  the  thickness  definition  is  not  required,  since  only  the  camber  surface  (de¬ 
fined  by  the  airfoil  mean  lines)  is  included  in  the  lifting  solution.  Since  the  flow  is 
supersonic,  a  simple  aft-marching  solution  can  be  used.  The  drag-due- to-lift  incre¬ 
ment  is  then  found  by  the  integration  of  the  pressure  force  components  on  the  camber 
surface  in  the  axial  direction,  and  the  lift  is  found  by  the  integration  of  the  component 
of  the  lifting  pressures  in  the  normal  direction. 

Attainable  vortex  (leading  edge)  thrust  coefficient  values  are  estimated  and 
included  in  the  drag-due-to-lift  calculation.  The  method  used  to  estimate  the 
attainable  vortex  thrust  was  developed  and  outlined  in  reference  8.  The  authors 
of  reference  8  used  .simple  sweep  theory  and  a  study  of  experimental  two-dimensional 
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airfoil  data  to  define  an  empirically  based  attainable  local  thrust  factor,  AY,  as  a 
function  of  the  airfoil  geometry,  Mach  number  normal  to  the  leading  edge,  and  the 
Reynolds  number.  AY  is  defined  as  the  ratio  of  the  attainable  leading  edge  thrust 
coefficient  to  the  value  corresponding  to  full  vortex  thrust  (fully  attached  flow).  The 
values  of  AY  vary  from  zero  for  vortex  flow  to  unity  for  fully  attached  flow.  The 
essential  ideas  of  this  theory  are  illustrated  in  Figure  2.4.  To  calculate  the  local 
attainable  thrust  factor,  AY,  for  a  given  wing  section  using  this  procedure  the  local 
leading  edge  radius  must  first  be  defined,  The  leading  edge  radius  of  a  given  airfoil  is 
determined  (reference  8)  by  fitting  a  square-root  function  to  the  thickness  distribution 
using  the  equation 

;  -T 

—  =  AY/-  (2.51 

2c  «  c 

where  the  constant  A’  can  be  evaluated  by  selection  of  some  f  near  the  leading  edge 
and  the  corresponding  Y  The  equation  for  a  circle  of  radius  ^  which  is  tangent  to 
the  half  thickness  curve  (or  the  airfoil  surface)  at  the  leading  edge  is 


(2.6) 


d  can  be  solved  with  equation  2.5  at  y=0  to  obtain  the  leading  edge  radius  ratio. 


r_=!0_ 

c  2 


(2.Y) 


Note  that  this  expression  gives  the  leading  edge  radius  in  the  chordwise  direction. 
Simple  sweep  theory  is  used  to  calculate  the  and  maximum  :  normal  to  the  leading 
edge,  which  are  used  with  the  Reynolds  number  and  Mach  number  normal  to  the 
leading  edge  to  calculate  AY  based  on  empirically  derived  equations. 


2.3  Friction  Drag 

Since  the  primary  purpose  of  HLFC  systems  is  to  reduce  the  friction  drag  ever  the 
wing,  this  is  a  critical  part  of  the  analysis.  To  get  a  single  correct  value  for  the  drag 
component  due  to  friction,  Cor  the  pressure  distribution  over  the  entire  wing  would 
be  found  using  an  Euler  code  (or  some  other  high-order  method).  This  pressure 
(or  velocity)  distribution  would  then  be  used  to  perform  a  boundary  layer  stability 
analysis  using  a  code  such  as  COSAL  (reference  9).  Knowing  what  types  of  flow 
exist  over  various  parts  of  the  wing,  the  corresponding  local  friction  coefficients  c/  for 
various  surface  regions  would  then  be  calculated  and  the  total  friction  drag  coefficient. 
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Cpn  of  the  wing  could  be  obtained.  An  even  more  accurate  but  very  time  consuming 
method  would  be  to  obtain  Xavier- Stokes  solutions  over  the  wing  which  would  include 
the  effects  of  viscosity. 

Performing  these  calculations  on  a  single  configuration  is  time  consuming,  and 
doing  so  for  a  large  number  of  configurations,  at  least  in  this  application,  is  prohibitive 
due  to  time  constraints,  Therefore,  some  assumptions  will  be  made  as  to  the  type  of 
flow  (laminar  or  turbulent)  that  should  east  over  the  surfaces  of  the  wings.  The  wings 
designed  by  the  camber  optimization  method  will  be  assumed  to  have  turbulent  flow 
over  the  entire  wing,  though  they  may  actually  be  capable  of  maintaining  small  areas 
of  laminar  flow  under  most  flight  conditions.  The  wings  designed  for  LFC  will  be 
evaluated  assuming  various  amounts  of  laminar  flow  over  their  upper  surfaces,  giving 
optimistic,  conservative,  and  most  likely  estimates  for  the  friction  drag. 

The  friction  drag  calculations  are  based  primarily  on  the  T'  method  for  the 
turbulent  flow  regions,  and  on  Blasius’  theory  for  the  areas  of  laminar  flow'.  Since 
both  of  these  methods  assume  flow  over  a  flat  plate,  the  components  of  the  aircraft  are 
represented  as  strips  of  appropriate  wetted  areas  and  reference  lengths.  Components 
like  the  wing,  which  vary  greatly  in  chord  length,  must  be  divided  spanwise  into 
several  strips  for  better  approximations  of  the  friction  drag.  The  drag  increments  on 
all  of  the  strips  are  then  summed  to  obtain  the  total  friction  drag  of  the  configuration. 


2.3.1  Turbulent  Boundary  Layer  Friction  Drag 

The  skin  friction  drag  for  turbulent  flow  is  computed  using  the  T'  method  outlined  by 
Sommer  and  Short  in  reference  10.  This  method  assumes  that  the  surface  is  a  smooth 
flat  plate  with  turbulent  boundary  layer  conditions.  Also,  transition  is  assumed  to 
occur  at  the  leading  edge  of  each  component  or  strip. 

The  T  method  is  based  on  the  calculation  of  a  compressible  skin  friction 
coefficient  C/,  from  a  reference  skin  friction  coefficient,  C'j,  for  the  given  Mach 
number,  Reynolds  number,  Rr and  adiabatic  wall  temperature.  T\yad.  The 

atmospheric  pressure,  p^,  temperature,  T, kinematic  viscocdty,  p,*.,  and  Reynolds 
number  per  Mach  per  unit  length,  are  first  calculated  for  the  given  flight  altitude 
using  1962  standard  day  atmosphere  equations. 

The  following  calculations  are  then  iterated  for  each  component  to  determine 
the  non-adiabatic  wall  temperature,  Tw,  and  skin  friction  coefficient,  Cj.  For  the  first 
iteration,  the  wall  is  assumed  to  be  adiabatic,  and  wall  temperature  is  approximated 
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using  the  compressible,  one-dimensional,  adiabatic  wall  equation 


—  1  -f  r 


-  1 


Ml 


(2.8) 


If  we  assume  that  7=1.4  for  air,  and  a  wall  recovery  factor  r  of  0.88,  we  get 

=  Id  0.176  Ml  (2.9) 

CO 

Sommer  and  Short  (reference  10)  obtained  the  key  relationship  to  this  theory  based 
on  experimental  data,  which  is 


~  =  1  *  0.035  Ml  ~  0.45  -  l)  (2.10) 

*  X  '  *  CO  ' 

The  Reynolds  number  relationship,  assuming  a  constant  pressure  through  the  bound¬ 
ary  layer  and  a  thermally  perfect  gas,  is 


R 


R 


OO 


(2.11) 


in  which  the  dynamic  viscosity  ratio  is  given  by  the  Sutherland  relationship  (with  the 
Sutherland  viscosity  constant  set  to  198.72  ’R) 


/  r\u  /Too  +  198. 72 \ 
[T^J  vr'+  198.72  / 


(2.12) 


where  the  T's  are  expressed  in  degrees  Rankine.  The  Reynolds  number  found  by 
equation  2.11  is  used  to  calculate  the  incompressible  skin  friction  coefficient,  given  by 
the  Karman-Schoenherr  equation 


logw{C/,Ri) 


which,  using  the  T1  analogy,  gives 


0.242 


logl0  (CjR) 


(2.13) 


(2.14) 


This  equation  must  be  solved  iteratively  to  obtain  Cj.  Then,  using  the  relationships 


Cf.  _  1 

C/  1  +  3.59 


(2.15) 
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and 


we  can  solve  for 


C  !_*  _ 

C,  ~  T 

j  i 


C, 


C)  To* 

1  +3.59 [ja,  T 


which  is  then  used  in  the  equation 


ft  a\ 

[it) 


(0.10425  x  10“12) 
n  C!t  (/?,o/0 


(2.16) 


(2-17) 


(2.18) 


Now'  the  non-adiabatic  wall  temperature  may  be  calculated  using 


y\v 


(2.19) 


Tliis  value  replaces  the  adiabatic  wall  temperature  as  an  initial  estimate  in  the  next 
iteration  which  starts  with  equation  2.10.  Once  equation  2.19  has  converged  (within 
five  iterations),  then  the  compressible  skin  friction  coefficient  C ;  is  found  using  the 
value  for  C's  from  equation  2.14  using  the  relationship 


C'j  Poo  V 


(2.20) 


2.3.2  Laminar  Boundary  Layer  Friction  Drag 


For  a  flat  plate,  if  viscous  effects  are  assumed  to  be  confined  to  a  thin  layer  over  the 
surface  and  a  uniform  velocity  is  assumed  outside  the  boundary  layer,  Blasius  found 
that  the  local  skin  friction  coefficient  at  a  distance  a;  behind  the  leading  edge  of  the 
plate  could  be  expressed  as 

0.664 

(2.21) 


Cl  = 


0.664 


where  the  local  Reynolds  number,  Rt  —  pVx/p.  Thus  the  average  skin  friction 
coefficient  Cj  for  one  side  of  the  flat  plate  of  unit  width  and  of  length  /  is  given  by 
the  relationship 


1.328 


(2.22) 


0 

where  /?;  —  pVl/p,  and  /  is  the  characteristic  length  of  the  plate  (reference  11). 
Excellent  agreement  between  theoretical  and  experimental  skin  friction  coefficients 
has  been  observed,  as  long  as  the  boundary  layer  actually  remains  laminar. 
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2.3.3  Mixed  Boundary  Layer  Friction  Drag 

As  previously  discussed,  the  T'  method  and  the  Blasius  method  were  used  together  to 
obtain  friction  drag  solutions  for  mixed  flows.  This  mixed  type  of  flow  is  assumed  to 
be  laminar  at  the  leading  edge  of  the  strip,  then  is  instantly  changed  to  fully  turbulent 
flow  at  a  specified  location  downstream  (Figure  2.5).  Of  course,  the  presence  of  the 
laminar  region  effects  a  lower  friction  drag  than  would  be  produced  by  a  fully  turbulent 
flow.  Following  the  teaching  of  L.  Prandtl,  this  reduction  in  drag  can  be  estimated  by 
assuming  that  behind  the  point  of  transition,  the  turbulent  boundary  layer  behaves 
as  if  it  were  turbulent  from  the  leading  edge.  Thus,  from  the  turbulent  friction  drag 
of  the  whole  length,  we  subtract  the  turbulent  friction  drag  of  the  part  up  to  the 
transition  and  add  the  laminar  friction  drag  for  that  same  length  (reference  12).  Ic 
is  crucial  that  these  drag  increments  be  based  on  the  correct  reference  lengths  in  all 
of  the  calculations.  The  drag  increment  produced  by  a  component  (strip),  defined  by 
a  wetted  area.  5u,ei,  a  reference  length,  L  and  a  specified  percentage  of  laminar  flow, 
x\ am,  can  be  calculated  using  the  relationship 

Co,  =  C,f„M  -  |cy,  -  cy,|i^  fat  (2.23) 

J-t]  l  ire/ 

where 


C/, 

CJ, 

i 

%cr\t 
3  wet 


turbulent  skin  friction  coefficient  of  the  entire  length  (based  on  /) 
turbulent  skin  friction  coefficient  of  the  laminar  part  (based  on  !„.,<) 
laminar  skin  friction  coefficient  of  the  laminar  part  (based  on  xcn <) 
the  reference  length  of  the  strip 

distance  from  the  leading  edge  at  which  transition  occurs 
wetted  area  of  the  strip 
reference  area 


The  turbulent  skin  friction  coefficients  C/,lur t  and  c/,  are  estimated  using  the 
T  method,  and  the  laminar  skin  friction  coefficient  C/,  is  calculated  using  the  Blasius 
formula.  The  total  skin  friction  drag  is  then  the  sum  of  the  component  drags 

.v 

CD,  =  YlCD,.n  (2-24 ) 

nr:  1 

where  ;V  is  the  toted  number  of  components. 
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A  computer  program  lias  been  written  to  calculate  the  Co,  for  mixed  flows 
based  on  these  combined  methods.  The  program  was  written  in  Fortran??,  and 
its  execution  times  were  very  sb  rt  due  to  the  quick  convergence  of  equations  2.14 
and  2.10.  The  wetted  area  calculations  were  made  with  a  code  which  previously 
existed  in  FortranIV.  but  was  converted  to  Fortran77.  This  code  breaks  down  a 
general  configuration  defined  in  the  Craidon  Wave-Drag  format  into  the  wetted  areas 
and  reference  lengths  required  by  the  skin  friction  program.  This  program  was  also 
modified  so  that  it  would  write  a  data  file  in  the  format  required  for  input  to  the 
skin  friction  program.  Flight  condition  cards  are  appended  to  the  end  of  this  file 
(as  many  as  desired)  which  specify  the  desired  Mach  numbers  and  altitudes  for  the 
calculations.  Using  this  file,  without  farther  modifications,  as  input  to  the  skin  friction 
<Lrag  program  results  in  the  calculation  of  the  configuration's  Co,  assuming  fully 
turbulent  flow  and  a  surface  emittance,  e,  of  0.8.  If  a  partially  (or  fully)  laminar 
solution  is  desired,  the  user  must  also  modify  the  input  file  to  include  the  desired 
amount  of  laminar  flow  on  the  upper  and  lower  surfaces  of  each  strip  of  the  wing. 
This  may  be  defined  by  either:  (1)  the  percent  area  (or  %  chord  length)  of  laminar 
flow;  (2)  the  specification  of  a  critical  Reynolds  number,  Br.criti  or  (3)  the  critical 
length,  icriti  at  which  the  flow  is  to  become  turbulent.  The  specified  parameter 
is  used  to  calculate  the  other  two  parameters  for  each  strip  at  each  Mach- altitude 
combination,  and  all  are  included  in  the  output.  The  program  output  also  includes 
the  average  skin  friction  coefficients,  C/,  and  skin  friction  drag  coefficients.  Con  at 
each  Mach-altitude  combination  for  each  component  (strip)  as  well  as  for  the  entire 
configuration. 


2.4  Validation  of  Pressure  Distribution  Estimates 

Since  the  programs  used  in  this  study  were  developed  elsewhere,  it  was  felt  necessary 
to  evaluate  the  validity  of  the  predictions  given,  and  the  present  application  of  the 
computer  programs.  The  theories  themselves  have  been  previously  evaluated  in 
references  13  and  14  and  have  since  been  incorporated  in  several  programs. 

The  Boeing  Design  and  Analysis  program  (to  be  discussed  in  chapter  3)  was 
used  to  estimate  all  pressure  distributions  presented  in  this  paper  and  is  evaluated 
here  for  accuracy.  Comparisons  between  theoretical  results  from  BDAP  and  the 
experimental  data  presented  in  reference  14  were  made  on  two  wings  of  the  same 
planfom  layout  (showu  in  Figure  2.6),  The  first  wing  was  not  cambered,  thus 
separating  the  effects  of  camber  in  the  solution.  A  second  wing  had  been  cambered 


11 


for  minimum  drag-due-to-lift,  using  optimization  techniques  which  were  available  at 
the  time,  at  a  design  lift  coefficient  of  0.08  (Figure  2.7). 

Evaluation  of  both  the  uncambered  and  cambered  wings  (Figure  2.8[a-g]) 
shows  the  theory  to  be  adequate  in  assessing  pressure  gradients,  and  even  pressure 
levels,  over  most  of  the  wing.  Agreement  between  the  predicted  and  experimental 
pressure  levels  appears  to  be  good  near  the  root,  but  it  deteriorates  near  the  wing 
tips,  especially  at  high  lift  conditions  (Figures  2.8[c]  and  2.8 [g] ) .  It  is  believed  that  this 
discrepancy  may  be  mostly  due  to  aeroelastic  effects  on  the  test  wing.  However,  while 
these  pressure  levels  near  the  tip  may  significantly  disagree,  the  pressure  gradients 
in  that  region  still  correspond  fairly  well.  The  predictions  of  the  pressure  levels  and 
gradients  near  the  leading  edge  also  prove  to  be  inferior,  especially  at  high  angles 
of  attack.  This  suggests  the  presence  of  a  vortex  flow  originating  at  the  wing  apex, 
producing  areas  of  separated  flow.  This  type  of  flow  is  not  predicted  by  the  linear 
methods  used;  thus,  its  effect  on  the  surface  pressure  is  not  observed.  These  pressure 
differences  introduce  little  error  in  the  prediction  of  the  aerodynamic  forces,  since 
the  areas  affected  are  small.  Indeed,  it  has  been  shown  in  references  13  and  14  and 
in  several  other  investigations  performed  since  then  (such  as  that  of  reference  6). 
that  when  the  force  estimates  are  combined  with  attainable  vortex  thrust  estimates 
(section  2.2),  the  force  predictions  correlate  very  well  with  test  data,  even  for  wings 
with  blunt  leading  edges.  The  small  errors  in  the  predicted  pressure  levels  can, 
however,  introduce  significant  errors  in  the  moment  predictions,  and  are  certainly 
a  concern  when  pressures,  as  in  the  present  case,  are  of  importance. 
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3  Designing  Wings  for 
Minimum  Drag-L  ue-to-Lift 


To  assess  the  effects  of  designing  wings  for  laminar  flow  control,  a  reference  design 
method  must  be  defined  upon  which  performance  comparisons  may  be  based.  The 
baseline  design  process  used  in  this  report  is  a  somewhat  standard  wing  loading 
optimization  scheme  that  is  presently  in  use  at  NASA's  Langley  Research  Center. 
Although  other  wing  optimization  methods  exist,  this  method  was  chosen  to  develop 
the  optimized  set  of  wings,  partially  for  convenience  since  it  is  used  by  the  Boeing 
Design  and  Analysis  (BDAP)  code  and  enabled  the  entire  parametric  analysis  to  be 
conducted  on  the  same  computer  system.  Examples  of  this  design  method  have  been 
illustrated  in  references  3-5.  The  technique  involves  the  determination  of  the  lifting 
load  distribution  which  will  produce  the  minimum  drag-due-to-lift,  and  the  definition 
of  the  camber  surface  required  to  produce  such  an  optimum  loading. 


3.1  Lifting  Load  Optimization  Theory 

The  wing  optimizations  in  this  paper  were  performed  using  the  wing  design  module  of 
the  Boeing  Design  and  Analysis  Program  (BDAP).  This  code  (references  15-18)  uses 
an  optimization  scheme  to  combine  several  basic  lifting  pressure  distributions  in  an 
optimum  fashion  in  order  to  obtain  a  total  wing  loading  distribution  which  produces 
the  minimum  drag-due- to-lift  at  the  design  lift  coefficient  (and  at  a  design  zero-lift 
moment  coefficient  if  desired).  The  lifting  pressure  distribution,  often  called  the  lift 
loading,  is  simply  the  difference  between  the  lower  and  upper  surface  pressures.  In 
linear  theory,  this  loading  is  assumed  to  be  entirely  independent  of  the  thickness 
distribution.  The  equations  for  the  candidate  lifting  loadings  which  are  considered 
for  the  design  of  the  wings  in  this  study  axe  listed  in  Table  3.1. 

The  performance  of  each  of  the  candidate  loadings  is  evaluated  by  integrating 
the  lifting  pressures  over  the  wing  to  obtain  its  overall  lift,  drag-due-to-lift,  and  zero- 
lift  pitching  moment.  In  addition,  the  interference  drag  increments  are  calculated. 
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Tabic  3.1:  Candidate  Loadings  Considered  in  Wing  Designs 


# 

&CP  = 

Loading  Description 

1 

k  a 

uniform  loading 

2 

hyb 

linear  spanwise 

3 

k  y 2 

quadratic  spanwise 

4 

k  iHn  - 1)1 

mid-span  loading 

5 

k  y/l-rj 

elliptical  spanwise 

6 

k  x'  c 

linear  chordwise 

7 

k  x1  2 

quadratic  chordwise 

8 

k  x'(c  -  x') 

parabolic  chordwise 

!  o 

k  x' 2 (1.5  c  —  x') 

cubic  chordwise 

10 

k  2(1  -r  5f )~° 5 

similar  to  flat  wing 

li 

&  ACptConj 

a  configuration  dependent 

12 

configuration  dependent 

at’s  =s  constants  of  proportionality. 
ky=spanwise  distance  from  wing  centerline. 
'Vrdongitudinal  distance  from  the  leading  edge. 


representing  the  drag  induced  on  each  loading  by  each  of  the  other  loadings.  Each 
loading  is  then  multiplied  by  a  constant  of  proportionality,  k,  such  fhat  its  resultant 
theoretical  lift  coefficient  is  unity.  This  eliminates  any  effect  that  machine  precision 
might  have  on  the  optimization  process  if  the  basic  loadings  were  to  produce  lift 
varying  greatly  in  magnitude.  These  aerodynamic  characteristics  are  supplied  to  the 
optimization  routine  where  Lagrange’s  method  of  undetermined  multipliers  is  used 
to  determine  the  weighting  factors,  M/1,  of  each  candidate  loading  that  will  minimize 
the  wing  axial  force  (drag-duc-to-lift)  while  producing  a  specified  normal  force  (lift) 
and,  if  requested,  a  specified  zero-lift  moment  coefficient.  During  the  optimization 
process,  these  aerodynamic  forces  axe  estimated  by 


and 


where 

n 

CL , 


Cl  =  EW<  Cl* 

i—  1 


cm  =  f2w< 

i=l 


(3.1) 


(3.2) 


number  of  candidate  loadings  fed  to  the  optimization  routine 
lift  coefficient  of  the  ith  loading 
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Cm.,  =  moment  coefficient  of  the  iih  loading 

Note  that  these  estimates  are  calculated  using  linear  theory  and  therefore 
include  no  attainable  vortex  thrust  considerations.  For  this  study  moment  constraints 
were  not  enforced  since  it  would  not  be  possible  to  apply  them  in  the  design  method 
of  Chapter  4,  During  preliminary  design  it  is  generally  considered  to  be  better  to 
trim  the  aircraft  by  movement  of  the  center  of  gravity  anyhow. 

Wing  loadings  which  are  independent  of  the  wing  itself  (such  as  body  buoyancy, 
body  upwash.  and  nacelle  buoyancy  pressures),  may  be  accounted  for  as  ' configuration 
dependent'  loadings  (candidate  loading  type  12  in  Table  3.1).  These  configuration 
dependent  loadings  are  forced  to  have  weighting  factors  of  unity  since  they  arc  not 
produced  by  the  wing  itself  and  thus  yield  no  camber  surface.  They  are  intended  to 
;dlow  the  optimization  of  the  camber  surface  in  the  presence  of  aerodynamic  wii  - 
loadings  caused  by  entities  other  than  the  wing  itself.  If  desired,  the  configurai,r 
dependent  loadings  may  also  be  used  as  candidate  loadings  in  the  optimization 
process  (loading  type  11)  to  provide  a  good  means  of  using  camber  to  capitalize 
on  any  beneficial  effects  of  an  external  loading  (type  12),  or  to  cancel  the  lift  of  the 
independent  effect,  if  that  be  the  optimum  solution.  In  such  an  optimization,  the 
configuration  dependent  loadings  may  be  calculated  by  the  program;  or,  they  may  be 
input  by  the  user  in  tabular  form. 

Since  linear  theory  contains  no  limitations  on  the  surface  pressure  coefficients, 
it  can  predict  physically  impossible  pressures  less  than  vacuum  pressure.  Experimen¬ 
tal  data,  however,  indicates  that  it  is  advisable  to  reject  theoretical  solutions  when 
predicted  suction  values  exceed  70-80%  of  vacuum  pressure  (reference  19).  To  help 
insure  a  physically  possible  design  lifting  pressure  loading,  upper  surface  pressure  con¬ 
straints  are  applied  during  the  wing  design  process.  The  surface  pressure  is  generally 
constrained  to  be  greater  than  ^  of  vacuum  pressure. or 


Cp  >  .7 


-Poo 


Qoo 


=  .7 


-1 


(3.3) 


To  design  with  this  constraint,  the  thickness  pressures  are  first  evaluated  using  the 
near-field  wave  drag  module  so  that  it  may  be  included  in  the  limiting.  The  program 
also  provides  options  of  adding  camber  constraints  and  pressure  gradient  constraints 
at  given  locations  on  the  wing. 


3.2  Definition  of  the  Camber  Surface 


Since  the  aerodynamic  performance  evaluations  provided  in  the  design  mode  are 
actually  the  linear  combinations  of  several  camber  surfaces,  small  numerical  errors 
that  are  introduced  for  each  surface  can  accumulate  into  a  larger  error.  Thus,  for 
an  analysis  of  the  final  design,  the  separate  loadings  are  consolidated  to  defile  the 
optimum  lifting  loading  definition,  and  the  lifting  pressure  at  any  point  on  the  wing 
is  taken  to  be 

n 

&Cp(x,y)  ='^T\V,  &Cp,t(x,y)  (3  4) 

1=  1 


Linearized  theory  is  then  used  to  calculate  the  wing  surface  shape  required  to 
support  the  optimized  lift  distribution  (reference  15)  by  the  integral  equation 


Oz 

dx 


ix,y)  =  ~~  ACp(x,y)  +  j  R(x  -  £,y  -  rj)  ACP(£,??)  d3rj  (3.5) 


where 

£  _  X~Z _  (3  6) 

0Hy  -  if  -  0J  -  -  iY 

where 

3  =  v A/2  —  1  (supersonic) 
s  =  x  variable  of  integration 
■1  —  y  variable  of  integration 

r  —  interval  of  integration  defined  by  the  surface  of  the 
wing  planform  within  the  Mach  cone  from  (x,t/). 


The  inversion  method  by  which  the  BDAP  uses  the  lifting  pressure  distribution 
to  define  the  corresponding  camber  surface  is  known  to  often  predict  camber  surfaces 
that  are  too  severe  near  the  plane  of  symmetry.  The  designer  may  either  avoid  this  by 
imposing  camber  constraints  at  the  root  of  the  wing  during  the  optimization  process, 
or  he  can  correct  it  by  substituting  a  milder  (or  even  flat)  camber  definition  at  the 
root  of  the  wing  subsequent  to  the  calculation  of  the  camber  surface.  This  camber 
surface  definition  may  then  be  combined  with  the  thickness  definition  and  evaluated 
using  whatever  methods  are  at  hand.  The  evaluation  techniques  used  in  this  paper 
have  already  been  discussed  in  detail  in  Chapter  2. 
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3.3  Application  to  a  Sample  Wing 


The  design  method  described  in  this  chapter  was  applied  to  a  baseline  planform  with 
A  =  ,1  and  AR^ 2.0  which  subsequently  will  be  referred  to  as  the  baseline  planform. 
The  planform  (reference)  area  of  the  wing  was  chosen  to  be  8000  square  feet  and  the 
leading  edge  sweep  was  chosen  to  be  67  *  .  The  wing  was  designed  for  a  Mach  number 
of  2.0  at  a  cruise  altitude  of  55,000  feet.  A  NACA-64  series  airfoil  thickness  profile 
with  a  maximum  thickness  ratio  of  4%  was  used  for  the  wing  (Figure  3.1).  The  basis 
for  these  selections  will  bo  discussed  in  Chapter  5. 

In  using  the  wing  design  module  of  BDAP.  several  inconsistencies  were 
observed.  After  evaluating  camber  surfaces  created  by  several  design  runs,  it  was 
found  that  using  all  of  the  candidate  loadings  (in  Table  3.1)  in  the  optimization  did 
not  necessarily  lead  to  the  best  performing  wing.  In  the  present  analysis,  use  of  the 
'oading  similar  to  a  flat  plate  (#  10)  was  found  to  cause  quite  poor  results.  It  was 
later  discovered  that  many  users  have  had  trouble  with  this  loading,  and  it  is  generally 
recommended  that  it  not  be  used.  This  may  be  related  to  the  fact  that  this  loading  is 
the  only  one  for  which  the  difference  between  the  leading  edge  A Cv  and  the  trailing 
edge  ACP  is  nonzero  and  in  no  way  a  function  of  the  local  chord  length.  It  can  be 
easily  found  that  for  this  loading,  lifting  pressure  at  the  leading  edge  is  (\/6)  times 
the  lifting  pressure  at  the  trailing  edge  for  any  value  of  c.  This  can  result  in  a  very 
large  gradient  in  A Cp  for  areas  where  the  local  chord  length  is  very  short,  causing 
liighly  tapered  wir.gs  to  be  excessively  twisted  at  the  tips. 

It  was  also  found  that  excluding  the  chordwise  cubic  loading  (i  9)  resulted 
in  wings  which  performed  almost  as  well  as  when  the  loading  was  included,  while 
producing  a  much  milder  camber  surface.  This  paradox  is  mentioned  in  reference  16, 
where  it  is  suggested  that  the  combination  of  the  higher  order  chordwise  loadings 
tend  to  produce  excessive  twist  or  waviness  in  the  calculated  camber  surface  unless 
constraints  are  imposed  on  the  wing’s  upper  surface  pressure  level  and  gradient.  It 
was  found,  through  several  design  calculations  involving  different  combinations  of 
candidate  loadings  and  constraints,  that  if  reasonable  pressure  gradient  constraints 
were  imposed  at  the  root  section  of  the  wing,  the  unwanted  waviness  in  the  calculated 
camber  surface  was  minimized,  even  when  using  loadings  5-8  simultaneously  in  the 
optimization.  An  overall  piessure  minimum  constraint  of  —  vacuum  pressure  was 
also  applied  during  the  design  process  but  was  found  to  be  unnecessary  due  to  the 
low  design  Cl,  of  the  wing,  Thus,  the  wing  was  designed  using  loadings  1-8  in  Table 
3.1  with  pressure  gradient  constraints  at  the  root. 
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The  most  heavily  weighted  chordwise  loadings  from  the  optimization  were 
found  to  be  the  linear  (#6)  and  quadratic  (#7)  candidates,  and  the  most  lightly 
weighted  was  the  parabolic  candidate  (#8).  It  was  also  observed  that  the  elliptical 
loading  (#5)  was  the  most  heavily  weighted  spanwise  loading,  and  the  mid-span 
loading  (#4)  was  the  most  lightly  weighted.  The  uniform  loading  (#1)  was  also 
among  the  most  heavily  weighted  loadings. 

The  baseline  planform.  designed  in  such  a  way,  is  shown  in  Figure  3.2  and 
its  corresponding  pressure  distribution  in  Figure  3.3.  As  previously  discussed  in  this 
chapter,  the  inverse  solution  in  the  BDAP  program  often  yields  excessive  camber  at 
the  root  of  the  wing  (as  seen  in  Figure  3.2).  This  was  corrected  here  by  linearly 
extrapolating  the  camber  values  from  the  next  two  span  stations  to  obtain  a  new 
root  camber  definition.  Softening  the  camber  actually  improved  the  performance 
of  the  wing  as  evaluated  by  the  methods  described  in  Chapter  2.  This  is  probably 
because  the  evaluation  methods  are  of  a  slightly  higher  order  than  the  design  methods, 
including  attainable  leading  edge  thrust.  The  resulting  wing,  for  our  baseline  planform 
is  shown  in  Figure  3.4,  and  its  airfoil  sections  are  shown  more  dearly  in  Figure  3.5.  As 
shown  in  Figure  3.6,  this  wing  design  possesses  a  much  flatter  upper  surface  pressure 
distribution  than  before  the  root  modification,  though  this  was  not  a  goal  in  the 
design  process. 


18 


4  An  Inverse  Design  Method 
for  HLFC  Wings 

4.1  Design  Goals 


The  objective  of  this  design  procedure  is  to  camber  the  wing  such  that  it  produces 
an  upper  surface  pressure  distribution  which  promotes  natural  laminar  flow  over 
the  mid- section  of  the  upper  surface,  thus  making  the  wing  amenable  to  a  hybrid 
laminar  flow  control  (HLFC)  system.  In  references  2D  and  21,  the  authors  discuss 
the  type  of  upper  surface  pressure  distribution  which  maximizes  the  design  Mach 
number  of  supercritical  laminar  flow  control  (SC  LFC)  airfoils  while  maintaining 
attached  flow.  A  pressure  distribution  with  a  low  supersonic  flat  rooftop,  preceded 
by  a  supersonic  front  pressure  minimum,  has  shown  itself  to  be  superior.  Since  linear 
potential  methods  are  to  be  employed  in  this  design  process,  the  pressure  rise  near 
the  trailing  edge  will  not  be  noticed.  Also,  since  this  method  was  shown  (in  Section 
2.4)  to  be  somewhat  deficient  in  estimating  the  pressures  near  the  leading  edge,  the 
small  front  pressure  minimum  will  not  be  sought.  With  these  approximations,  the 
pressure  distribution  that  we  seek  on  the  upper  surface  simplifies  to  a  constant  in  the 
chordwise  direction,  often  called  a  'fiat-roof-top  pressure  distribution. 

To  design  the  camber  surfaces  of  the  laminar  flow  wings,  the  Boeing  Design  and 
Analysis  Program  (BDAP)  described  in  Chapter  3  is  incorporated.  This  program  is 
not  intended  to  directly  design  a  wing  for  any  given  upper  surface  pressure  distribution 
since  it  bases  the  camber  surface  designs  on  the  lifting  pressure  distribution.  This 
necessitates  a  somewhat  unconventional  use  of  the  code,  but  as  will  be  shown,  it  allows 
the  design  process  to  include  an  optimization  of  the  spanwise  lifting  distribution  for 
minimum  drag-due-to-lift. 
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4.2  Inverse  Design  Procedure 


This  method  takes  advantage  of  the  linear  theory  assumption  that  if  a  camber  surface 
can  be  expressed  as  a  linear  combination  of  two  or  more  camber  portions,  then  the 
aerodynamic  characteristics  of  the  camber  surface  may  evaluated  as  the  sura  of  the 
aerodynamic  characteristics  of  these  camber  portions.  The  lifting  loading  of  the  HLFC 
wing  will  be  considered  to  be  made  up  of  two  portions  which  will  be  linearly  combined 
to  define  the  design  wing  loading.  The  first  portion  will  be  referred  to  as  the  flat- 
roof-top  loading.  This  loading  is  one  which  will  combine  with  the  thickness  pressures 
to  produce  a  flat  chordwise  upper  surface  pressure  distribution,  while  producing  little 
or  no  lift  across  the  span.  In  other  words,  it  is  desired  that  the  section  lift  coefficient 
at  each  span  station  and  the  total  lift  both  be  zero.  This  type  of  lifting  pressure 
distribution  is  illustrated  in  Figure  4.1  for  a  wing,  and  in  Figure  4. 3 [a]  for  a  given 
span  station  on  the  wing. 

The  second  portion  of  the  total  lifting  loading,  which  is  responsible  for 
producing  the  wing  lift,  will  be  called  tht  3pan- optimized  loading.  This  name  was 
chosen  since  this  loading  is  to  have  a  spanwise  lifting  distribution  which  is  optimized 
for  minimum  drag-due-to-lift,  while  having  a  constant  chordwise  lifting  pressure  at 
each  span  station  (i.e.  pure  aerodynamic  twist).  The  optimization  is  performed  ir. 
the  presence  of  the  flat-roof-top  loading.  The  span-optimized  loading  is  illustrated  in 
Figure  4.2  for  a  wing,  and  in  Figure  4.3(b]  for  a  given  span  station  on  the  wing. 

These  two  pressure  distributions  will  combine  with  the  thickness  distribution 
(Figure  4.3(c))  to  define  a  wing  surface  pressure  distribution  which  possesses  a  flat- 
roof-top  in  the  chordwise  direction  and  a  spanwise  lifting  distribution  optimized 
for  minimum  drag-due-to-lift  (Figure  4.3(d)).  The  camber  surface  corresponding 
to  each  of  these  loadings  is  calculated  using  linear  theory.  Once  these  camber 
surface  descriptions  axe  defined,  they  are  combined  with  the  thickness  description 
and  evaluated  to  see  how  well  the  total  upper  surface  pressure  gradients  correspond 
with  those  desired.  If  these  are  satisfactory,  the  wing  is  ready  to  be  evaluated  for  its 
overall  aerodynamic  characteristics  using  the  methods  described  in  Chapter  2. 


4.2.1  Determination  of  the  Flat-Roof- Top  Loading  Distribu¬ 
tion 

Since  the  BDAP  code  uses  lifting  pressures  (loadings)  as  a  basis  for  the  optimization 
of  the  design,  it  is  necessary  to  determine  the  loading  distribution  which  would 


combine  with  the  thicKness  pressures  to  produce  the  desired  upper  surface  pressure 
distribution,  This  loading  distribution  is  dependent  on  the  planform  shape  as  well  as 
the  thickness  profile  of  the  wing, 

The  lifting  pressure  coefficient  is  defined  as  the  difference  between  the  lower 
and  upper  surface  pressures,  or: 


ACP  —  CPl  CPu 

(4.1) 

In  linear  theory,  the  upper  and  lower  surface  pressures 
combination  of  the  effects  of  tliickness  and  camber. or 

are  assumed  to  be  a  linear 

-  (?*.,  +  Cp„.£ 

(4.2) 

mid 

Cpt  =  Cpi,<  +  CPl  c 

(4.3) 

where 


CPu  =  lower  surface  pressure  coefficient  due  to  camber 
CPu  c  =  upper  surface  pressure  coefficient  due  to  camber 
CP| ,  -  lower  surface  pressure  coefficient  due  to  thickness 
CPu ,  =  upper  surface  pressure  coefficient  due  to  thickness 


Due  also  to  the  linearity  of  the  solution,  the  pressures  due  to  camber  follow  the 
relationship 

Cp..c  =  ~Cpu 

and  the  pressure  due  to  thickness  is  the  same  on  the  top  surface  as  on  the  bottom,  or 


so 


Cpt  —  Cp„ ,  —  CPl , 

(4.5) 

+  Cp,.c  —  CP<1 ,  —  CPu  e  =  —2  CPu  c 

(4.6) 

gn  case  : 

Pu.de*  ^Pt i,l  ”b  ^-'Pu.c.de* 

(4.7) 

From  linear  theory,  the  surface  pressure  due  to  thickness  on  the  wing  is  simply  the 
average  of  the  upper  and  lower  surface  pressures,  or 

Cpu  +  Cpi 


CPM  — 


(4,8) 
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so 


C, 


£pu  +  Cp\ 


Pu.c.dei 


2  CpuMtl 


and  using  equation  4.6 


A r  -  1  l'Cpu  ~  Cp>-  -  r 

^^Pde,  ~  1  {  o  ^Pu.dtt 


(4.9) 


(410) 


For  a  given  wing  planform,  with  a  given  thickness  distribution,  the  upper 
and  lower  surface  pressures  (CPu  and  CPl )  are  calculated  at  any  lift  condition  for  the 
uncambered  wing.  These  pressures  are  then  used  in  Equation  4.10  to  determine  the 
loading  distribution  (A CPdti)  required  to  sustain  any  given  pressure  distribution  over 
the  upper  surf.’  ce  of  the  wing  {CPudtl). 

For  our  case,  we  are  not  really  interested  in  maintaining  a  specific  upper  surface 
pressure  distribution,  but  rather  a  specific  upper  surface  pressure  profile,  or  gradient 
distribution.  The  profile  we  seek  has  been  simplified  to  be  a  flat  one  (flat-roof¬ 
top).  Thus,  any  chordwise  constant  CPu  dti  may  be  used  in  equation  4.10  to  give 
an  equally  suitable  upper  surface  pressure  distribution.  To  design  the  flat-roof- top 
camber  portion,  a  constant  Cpu,4„  is  found  at  each  span  station  such  that  the  average 
lifting  pressure  loading  (ACP)  is  nearly  zero,  or 

l 

J  d(-)  ~  0  (411) 

Q 

This  is  done  in  an  attempt  to  maintain  a  spanwise  loading  distribution  and  overall 
lift  coefficient  of  nearly  zero  for  this  camber  portion.  A  computerized  spreadsheet 
was  found  to  be  very  useful  in  finding  the  CPudtt  required  to  do  this,  and  in  creating 
the  input  data  records  for  the  wing  design  program. 

4.2.2  Determination  of  the  Span-  Optimized  Loading  and 
Camber  Surface 

In  the  design  of  this  surface  portion,  the  optimization  capabilities  of  the  BDAP  code 
are  utilized  The  flat-roof-top  loading  distribution  is  included  in  the  optimization  as  a 
configuration  dependent  loading  (loading  type  12  in  Table  3.1).  The  BDAP  program 
is  used  to  determine  the  optimum  combination  of  the  candidate  lifting  loadings  which 
were  constant  in  the  chordwise  direction  (loadings  1-5  in  Table  3.1)  in  the  presence 
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of  the  flat-roof-top  loading.  Since  the  flat-roof-top  loading  is  not  considered  by  the 
program  to  be  a  camber-dependent  loading,  its  weighting  factor  is  forced  to  be  equal  to 
unity  in  the  optimization  process  and  is  not  included  in  the  camber  surface  calculation. 
Linear  theory  is  then  used  to  calculate  the  span-optimized  portion  of  the  camber 
surface  from  the  span-optimized  portion  of  the  loading  distribution.  When  the  span- 
optimized  portion  of  the  lifting  surface  is  evaluated  alone,  its  loading  distribution  is 
seen  to  be  constant  in  the  chordwise  direction  and  somewhat  elliptical  in  the  spanwise 
direction. 


4.2.3  Determination  of  the  Flat*Roof-Top  Camber  Surface 

The  portion  of  the  camber  surface  required  to  produce  the  flat-roof-top  loading 
distribution  is  obtained  using  BDAP  by  bypassing  its  optimization  routine.  The 
program  is  used  to  find  the  weighting  factor,  IV,,  which  is  required  to  enable  a  single 
candidate  loading,  proportional  to  the  flat-roof-top  loading  (input  as  loading  11  in 
Table  3.1),  to  meet  a  specified  lift  coefficient  criterion.  When  such  a  non-optimized 
solution  is  requested,  the  program  will  calculate  the  camber  surface  required  to 
support  even  a  configuration-dependent  loading. 

To  know  what  the  design  lift  coefficient  should  be  in  this  computation, 
the  results  obtained  during  the  calculation  of  the  span-optimized  loading  must  be 
used.  During  this  process  each  candidate  loading  was  evaluated  for  its  aerodynamic 
properties  before  being  submitted  to  the  optimization  routine.  This  was  also  true 
for  the  flat-roof-top  loading,  since  its  lift  increment,  denoted  by  Ci}rt  ,  had  to  be 
accounted  for  in  the  total  lift.  When  Cl,  r  t  is  used  as  our  design  lift  coefficient 
and  the  optimization  routine  is  bypassed,  the  weighting  factor,  IV,,  required  to  meet 
this  criterion  should  be  near  unity  since  this  loading  alone  is  expected  to  meet  the 
design  lift  condition.  The  camber  definition  required  to  support  the  input  loading 
distribution  times  this  weighting  factor  is  then  calculated  by  the  program  using 
linearized  theory  (equation  3.5). 


4.2.4  Total  Camber  Surface 

Using  the  linear  theory  assumptions  made  at  the  beginning  of  this  chapter  the 
flat-roof-top  portion  and  the  spanwise-optimized  portion  of  the  camber  surface,  are 
combined  to  define  the  design  camber  surface.  This  camber  surface  will  produce 
a  spanwise  lift  distribution  designed  for  minimum  drag-due-to-lift,  and  produce  a 
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cliordwise  lifting  pressure  distribution  which  will  combine  with  the  thickness  pressures 
to  yield  a  flat-roof-top  pressure  distribution.  As  discussed  in  Chapter  3,  the  linear 
theory  inversion  method  used  by  BDAP  calculates  camber  surfaces  that  are  usually 
too  severe  near  the  plane  of  symmetry.  In  a  typical  optimization  design  process,  the 
designer  has  two  options  to  correct  this.  He  may  either  impose  camber  constraints  at 
the  root  of  the  wing  during  the  optimization  process,  or  he  can  substitute  a  milder 
camber  definition  at  the  root  of  the  wing  subsequent  to  the  calculation  of  the  camber 
surface,  Since  our  solution  to  the  flat-roof-top  portion  of  the  design  already  has  one 
constraint  (lift  coefficient)  and  only  one  loading  distribution  to  satisfy  it  with,  no 
further  constraints  are  possible.  For  this  reason,  after  the  design  camber  surface  is 
obtained  using  the  methods  previously  outlined,  the  root  section  should  be  replaced 
with  a  much  milder  definition.  It  was  felt  that,  for  this  use  of  the  program,  it  would 
be  best  to  rotate  the  root  airfoil  until  the  three  inboard  trailing  edge  camber  values 
formed  a  straight  line.  That  is,  the  three  inboard  sections  would  have  a  linear 
relationship  in  twist.  In  this  way,  the  calculated  camber  definition  of  the  inboard 
airfoil  section  would  be  maintained. 


4.2.5  Summary  of  Inverse  Design  Procedure  for  HLFC 
Wings 

(a)  Define  the  wing  planform  geometry  and  thickness  distributions. 

(b)  Evaluate  the  flat  wing  for  pressure  distribution. 

(c)  Use  the  flat  wing  pressures  in  the  equation 

AC,..  =  2  (Z±±ClL  _  c^..)  (4-10) 

adjusting  CPu at  each  span  station  such  that 

l 

J  A *0  (4.11) 

o 

to  define  the  flat-roof-top  loading  distribution. 

(d)  Calculate  the  span-optimized  loading  distribution,  using  only  the  chordwise- 
constant  candidate  loadings  in  BDAP’s  wing  design  module,  in  the  presence 
of  the  flat-roof-top  loading.  This  procedure  also  yields  the  calculated  lift 
increment  produced  by  the  flat-roof-top  loading,  Cijrt  . 
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(e)  Determine  the  camber  surface  required  to  produce  the  flat-roof-top  loading 
using  BDAP’s  wing  design  module.  This  is  done  by  using  the  flat-roof-top 
loading  as  the  only  candidate  loading,  bypassing  the  optimization  routines, 
and  designing  at  Citr,  . 

(f)  Combine  the  flat-roof-top  and  span-optimized  camber  surface  definitions  to 
define  the  design  HLFC  camber  surface. 

(g)  Adjust  the  twist  definition  at  the  root  to  provide  a  more  continuous  spanwise 
camber  definition. 

4.3  Application  to  a  Sample  Wing 

The  design  procedure  described  in  this  chapter  was  used  to  develop  an  HLFC  wing 
having  the  same  planform  and  thickness  distribution  as  the  wing  which  was  designed 
for  minimum  diag-due-to-lift  using  the  method  of  Chapter  3. 

As  discussed  previously,  this  design  procedure  involves  the  definition  of  the 
design  camber  surface  in  two  portions.  The  uncambered  wing  was  first  evaluated 
to  determine  its  upper  and  lower  surface  pressure  distribution.  This  distribution  is 
shown  in  Figure  4  4  at  the  design  lift  coefficient  of  0.08  for  the  baseline  planform.  This 
pressure  distribution  was  then  used  in  equation  4.10  to  find  the  flat -roof-top  loading 
distribution.  This  is  the  loading  that  combines  with  the  thickness  pressures  to  create 
a  constant  chordwise  upper  surface  pressure  distribution.  The  design  upper  surface 
pressure,  CPiWe>,  in  equation  4.10  was  specified  to  be  a  constant  value  at  each  span 
station  and  was  varied  until  the  average  A Cp  was  nearly  zero,  thus  approximating 
the  zero  lift  condition. 

The  wing  design  module  of  the  BDAP  program  was  then  used  to  calculate 
the  span-optimized  loading  distribution  and  the  corresponding  camber  surface.  This 
was  achieved  by  optimizing  the  candidate  loading  distributions  which  are  constant 
in  the  chordwise  direction  (loadings  1-5)  in  the  presence  of  the  flat-roof-top  loading. 
As  was  the  case  for  the  wing  designed  for  minimum-drag-due-to-lift  (Section  3.3),  it 
was  observed  that  the  elliptical  loading  (#5)  was  the  most  heavily  weighted,  and  the 
mid-span  loading  (#4)  was  the  most  lightly  weighted. 

Next,  the  flat-roof-top  camber  surface  was  determined  by  bypassing  the 
optimization  routine  of  BDAP  and  simply  finding  the  surface  required  to  obtain 
the  flat-roof-top  loading  distribution.  In  this  procedure,  the  lift  coefficient  increment 
produced  by  the  flat-roof-top  loading,  Cl,  r  (  ,  was  used  as  the  design  lift  coefficient. 
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Table  4.1:  Lift  Contributions  of  Surface  Portions  for  Wings  in  Parametric  Study 


Wing 

ID 

A 

i 

AR 

CLfr, 

flat-roof-top 

portion 

i 

^ Ltpan  —  opl 

span-optimized  i 
portion 

W015 

0.0 

1.5 

-0.00836 

0,08836 

W020 

0.0 

2.0 

-0.00756 

0.08756 

|  W025 

0.0 

2.5 

-0.00746 

0.08746 

|  W115 

0.1 

1.5 

-0.00831 

0.08831 

W120 

0.1 

2.0 

-0.00761 

0.08761 

W125 

0.1 

2.5 

-0.00765 

0,08765 

W215 

0.2 

1.5 

-0.00841 

0.08841 

W220 

0.2 

2.0 

-0.00777 

0.08777 

W225 

0.2 

2.5 

-0.00777 

0.08777 

W315 

0.3 

1.5 

-0.00858 

0.08858 

W320 

0.3 

2.0 

-0.00812 

0.08812 

W325 

0.3 

2.5 

-0.00816 

0.08816 

Ci, , ,  had  been  identified  by  the  design  program  in  the  spanwise-optimization  process 
and  should  be  near  zero. 

The  pressure  distributions  exhibited  by  the  two  wing  portions,  including  the 
thickness  pressures,  are  shown  in  Figures  4.5  and  4.6  for  the  baseline  planform.  In 
these  figures,  the  two  are  evaluated  at  their  respective  lift  coefficient  increments, 
which  are  defined  by  the  relationship 

CLit,  =  CL 

tpan  —  opt  +  cL/rt. 

Thus,  the  sum  of  the  lift  coefficients  in  Figures  4.5  and  4.6  is  equal  to  the  design  lift 
coefficient  0.08.  The  lift  increments  of  the  two  camber  surface  portions  for  all  of  the 
wings  presented  in  this  paper  are  listed  in  Table  4.1. 

Once  the  two  portions  of  the  camber  surface  were  determined,  they  were 
combined  linearly  with  the  thickness  distribution  to  define  the  final  flat-roof-top 
spanwise-optimized  wing  (Figure  4.7).  This  wing  was  then  evaluated  and  the  total 
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pressure  distribution  is  shown  in  Figure  4,8. 

As  previously  suggested,  the  airfoil  section  at  the  root  of  each  wing  designed 
by  this  procedure  was  rotated  such  that  the  the  three  most  inboard  sections  would 
have  a  linear  relationship  in  twist.  The  baseline  HLFC  wing,  after  being  modified  in 
such  a  way,  is  shown  in  Figure  4.9,  and  the  airfoil  sections  are  illustrated  in  Figure 
4.10.  The  wing  with  the  modified  twist  at  the  root  was  then  re-evaluated,  and  as 
shown  in  Figure  4.11,  the  pressure  distribution  was  not  significantly  degraded  by  this 
modification. 
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5  Parametric  Study  of  the 
Wing  Design  Methods 


A  parametric  analysis  involving  several  arrow-wing  planforms  was  conducted  in  an 
attempt  to  identify  the  effects  that  aspect  ratio,  AR .  and  taper  ratio.  A.  may  have 
on  the  desirability  of  designing  wings  for  HLFC.  The  study  was  conducted  at  a  Mach 
number  of  2.0,  since  the  effect  of  camber  on  minimum  drag-due- to-lift  decreases  at 
lugher  Mach  numbers.  A  family  of  twelve  wings  of  varying  planform  shape  was 
developed  with  a  common  refeience  area,  leading  edge  sweep,  and  airfoil  thickness 
profile.  As  in  the  baseline  example  of  Chapters  4  and  5,  the  reference  planform  area 
was  chosen  to  be  8000  ft2  for  all  of  the  wings.  The  leading  edge  sweep  was  chosen  to 
be  67  * ,  placing  the  entire  wing  at  least  7  °  behind  the  Mach  cone.  This  results  in  a 
subsonic  Mach  number  of  0.78  normal  to  the  leading  edge.  Ln  this  study,  the  same 
NACA-64  series  airfoil  thickness  profile  with  a  maximum  thickness  ratio  of  4%  that 
was  employed  in  the  baseline  design  was  used  on  all  of  the  wings.  Arpec(  ratio  was 
varied  from  1.5  to  2.5  and  taper  ratio  was  varied  from  0.0  to  0.3.  This  variation  was 
selected  to  envelop  the  range  of  simple  arrow-wing  shapes  of  interest  for  supersonic 
HLFC.  The  resulting  family  of  planforms  is  shown  in  Figure  5.1  and  their  dimensions 
are  listed  m  Table  5.1. 

Each  of  these  planforms  was  cambered,  as  the  baseline  wing  was.  using 
both  the  minimum  drag-due-to-lift  method  described  in  Chapter  3,  and  the  HLFC 
design  procedure  of  Chapter  4.  The  aerodynamic  performances  of  the  wings  were 
then  calculated  using  the  methods  outlined  in  Chapter  2,  and  the  results  of  these 
evaluations  follow. 

5.1  Zero  Lift  Wave  Drag 

The  wave  drag  of  each  wing  was,  as  discussed  in  Chapter  2,  assumed  to  be  independent 
of  the  camber  definition.  Thus,  the  wave  drag  of  each  cambered  wing  could  be 
approximated  by  the  value  for  an  uncambered  wing  of  the  same  thickness  distribution, 
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Table  5.1:  Dimensions  of  Planforms  in  Parametric  Study 

a 

j  Wing  '  ;  T  | 


ID 

A 

All 

b 

Cr 

c, 

5 

Vi 

xi  ! 

W015 

0.0 

1.5 

109.54 

146.06 

0.00 

97.37 

18.26 

67.35 

W020 

0.0 

2.0 

126.49 

126.49 

0.00 

84.33 

21.08 

W025 

0.0 

2.5 

141.42 

113.14 

0.00 

75.42 

23.57 

74.38 

W115 

0.1 

1.5 

109.54 

132.78 

13.28 

89.33 

19.92 

69.25 

W120 

0.1 

2.0 

126.49 

114.99 

11.50 

77.36 

23.00 

73.52 

W125 

0.1 

2.5 

141.42 

102.85 

10.29 

69.19 

25.71 

77.87 

W215 

0.2 

1.5 

109.54 

121.72 

24.34 

83.85 

21.30 

71.14 

W220 

0.2 

2.0 

126.49 

105.41 

21.08 

72.62 

24.60 

VV225 

0.2 

2.5 

141.42 

94.28 

18.85 

64.95 

81.02 

W315 

0.3 

1.5 

109.54 

112.35 

33.71 

80.09 

22.47 

72.96 

W320 

0.3 

2.0 

126.49 

29.19 

69.36 

25.95 

78.47 

W325 

0.3 

2.5 

141.42 

26.11 

83.85 

°AU  dimensions  in  feet. 

Common  Dimensions  for  All  Wings: 

S  =  8000  fi\  Ait  =  67°,  =  <%,  r~  =  .1136%. 
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allowing  for  the  wave  drag  to  be  calculated  for  only  twelve  'vings  (shown  in  Figure 
5.1). 

Analysis  showed  that  the  wave  drag  decreased  with  increases  of  both  taper 
ratio  and  aspect  ratio  (Figure  5.2).  This  is  attributed  mostly  to  the  reduced  volume 
of  the  wings  with  higher  \  and  AR.  The  reduction  in  near-field  wave  drag  associated 
with  increasing  A  from  0.0  to  0.3  varied  from  7-14%  at  AR\  of  1.5  and  2.5,  respectively. 
The  reduction  in  wave  drag  associated  with  increasing  AR  from  1.5  to  2.5  also  varied 
from  7-14%  at  A’s  of  0.0  and  0.3,  respectively.  The  maximum  difference  in  wave 
drag  seen  between  any  two  wings  was  about  20%.  .All  of  the  values  that  have  been 
quoted  here  were  obtained  from  the  near-field  analysis,  but  as  shown  in  Figure  5.3. 
the  far-field  method  provided  very  similar  results. 


5.2  Skin  Friction  Drag 

The  turbulent  skin  friction  drag  was  also  assumed  to  have  no  direct  dependence  on 
the  camber  definition,  thus  the  small  effect  that  camber  has  on  the  wetted  area  of  the 
wing  is  neglected.  However,  for  wings  which  are  not  assumed  to  have  frilly  turbulent 
boundary  layers  (i.e.  those  designed  for  HLFC),  the  friction  drag  is  very  dependent 
on  the  camber  definition  by  means  of  its  influence  over  the  boundary  layer  tv., 
(laminar,  turbulent,  or  mixed).  In  this  study,  the  wings  which  had  been  optimi. 
for  minimum  drag-due- to-lift  with  no  boundary  layer  considerations  were  assumed 
to  possess  fully  turbulent  boundary  layers.  The  wings  that  were  designed  by  the 
HLFC  design  method  to  produce  flat- roof- top  pressure  distributions  were  evaluated 
assuming  varying  amounts  of  laminar  flow  on  the  upper  surface  and  fully  turbulent 
flow  on  the  lower  surface. 

On  a  simple  HLFC  wing,  suction  panels  could  extend  along  the  entire  span 
of  the  wing  for  some  given  ranges  in  percent  chord  length.  For  example,  an  HLFC 
wing  might  incorporate  one  suction  panel  from  the  leading  edge  to  10%  chord  length 
and  then  another  panel  between  30%  and  40%  chord  length.  The  flow  would  then 
be  stabilized  for  the  fust  10%  by  suction  and  a  region  of  natural  laminar  flow  would 
follow  until  the  flow  was  again  stabilized  by  the  second  suction  strip  at  30%  of  the 
local  chord  length.  Another  region  of  natural  laminar  flow  would  then  extend  for  some 
length  behind  the  second  strip.  If  it  is  assumed  that  the  local  mass  flow  rate  through 
a  suction  panel  is  proportional  to  its  width,  and  that  the  degree  of  stabilization 
accomplished  is  proportional  to  this  mass  flow,  then  the  increase  of  the  local  critical 
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Reynolds  number  due  to  the  suction  could  be  assumed  to  be  proportioned  to  the 
local  panel  width.  On  an  HLFC  wing  such  as  the  one  described  above,  the  local 
width  of  each  suction  panel  is  proportional  to  the  local  chord  length,  so  the  degree 
of  stabilization,  or  length  of  natural  laminar  flow  following  a  given  suction  panel 
is  proportional  to  the  local  chord  length.  On  this  basis,  a  study  was  conducted  in 
which  the  fraction  of  the  chord  length  which  was  laminar  was  assumed  to  be  constant 
over  the  entire  span  of  the  HLFC  wings.  The  results  of  this  study  are  discussed  in 
Section  5.2.1.  Calculations  were  also  made  assuming  a  constant  upper  surface  critical 
Reynolds  number  over  the  span  of  the  vring.  which  will  be  discussed  in  Section  5.2.2. 

5.2.1  Co,  Based  on  Critical  ~ 

Friction  drag  estimates  for  the  baseline  wing  (as  calculated  by  the  method  of  Section 
2.3)  are  shown  in  Figure  5.4  assuming  several  quantities  of  laminar  flow  on  the  upper 
surface  in  percent  chord  length,  where  the  skin  friction  drag  is  plotted  as  a  function 
of  Mach  number.  The  results  indicate  that  at  this  flight  condition,  if  laminar  flow 
could  be  maintained  over  the  entire  upper  surface,  there  would  be  approximately  a 
45%  reduction  in  skin  friction  drag  as  compared  to  a  fully  turbulent  wing.  It  can  also 
be  seen  that  the  relationship  between  the  amount  ot  Aminat  flow  and  the  friction 
drag  is  nearly  linear  lor  any  given  Mach  number.  This  is  also  true  for  the  other  wings 
in  'he  study. 

Figures  5.5  and  5.6  show  the  skin  friction  drag  at  cruise  for  Loth  fully  turbulent 
vi  and  wings  having  upper  surfaces  which  are  entirely  laminar.  The  figures  show 
that  for  the  fully  laminar  and  fully  turbulent  wings  friction  drag  increases  with  taper 
ratio  and  also  with  aspect  ratio.  The  increase  with  taper  ratio  appears  to  be  relatively 
small,  having  a  maximum  effect  on  friction  drag  of  less  than  2.5%  for  any  given  aspect 
ratio.  Increasing  the  aspect  ratio  from  1.5  to  2.5  has  slightly  more  influence,  increasing 
the  fricuon  drag  by  only  3. 5-4. 2%.  A  more  prominent  feature  of  this  figure  is  that 
the  difference  in  the  drag  coefficient  values  caused  by  a  change  in  AR  is  constant  for 
all  of  the  wings  with  a  given  A.  The  difference  between  the  AR  =  1.5  and  AR  —  2.5 
curves  is  about  .00011  for  the  turbulent  wings,  and  .00007  for  the  wings  with  fully 
laminar  upper  surfaces. 

Laminar  flow  can  probably  be  induced  over  no  more  than  75%  of  the  chord 
length  on  the  upper  surface  without  requiring  a  very  high  mass  flow  through  the 
suction  panels.  In  any  case,  beyond  this  point  several  physical  disturbances  such 
as  flap  hinge  Hues  would  probably  exist  on  an  HST  and  would  surely  trip  the 
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boundary  layer  unless  further  stabilization  methods  were  employed.  This  assumed 
value  corresponds  to  a  friction  drag  reduction  of  about  34%.  However,  it  will  also  be 
assumed  that  a  7°  turbulent  wedge  exists  at  the  wing-fuselage  abutment.  This  value 
is  based  on  extrapolated  subsonic  data,  as  jpven  by  Fisher  in  reference  22.  If  there  were 
any  wing  leading-edge  discontinuities,  the  area  bounded  by  a  14  °  wedge  with  its  apex 
at  the  discontinuity  would  also  be  assumed  turbulent  due  to  the  span  wise  disturbance. 
This  assumption  was  applied  to  all  of  the  wings,  and  new  friction  drag  estimates  were 
■.btained.  A  wing  with  a  higher  A  will  have  less  of  its  surface  area  submerged  in  the 
turbulent  wedge  than  a  wing  having  a  low  A  (holding  aspect  ratio  constant).  The 
introduction  of  the  wing  root  turbulent  wedge  into  the  calculations  greatly  reduced 
the  influence  of  taper  ratio  on  friction  drag  (Figure  5.7).  The  reduction  in  turbulent 
wedge  area  is  even  more  pronounced  with  an  increase  in  AR.  The  uppei  series  of 
curves  m  Figure  5.7  represent  the  wings  w.  h  75%  laminar  flow  on  the  upper  surface 
and  a  7  0  turbulent  wedge  originating  at  the  apex.  The  figure  shows  that  the  turbulent 
wedge  assumption  actually  reverses  the  effect  of  A  and  reduces  the  combined  effects 
of  A  and  AR  on  skin  friction  drag  to  less  than  2%.  For  the  planforms  with  low  As, 
the  influence  of  AR  is  essentially  eliminated  altogether.  When  this  case  is  compared 
to  the  fully  turbulent  wings,  friction  drag  reductions  of  about  28-32%  are  observed. 


5.2.2  Cdj  Based  on  Critical  Re 

If  it  were  assumed  that  the  flow  were  to  transition  from  laminar  to  turbulent  at  some 
specified  critical  Reynolds  number  as  is  assumed  on  many  natural  LFC  applications, 
mstead  of  at  a  given  percent  it  would  surely  change  the  predicted  dependence 
of  friction  drag  on  AR.  A  fven  Rta-n  would  translate  into  a  much  smaller  area  of 
laminar  flow  on  a  low  AR  (short  span)  wing  than  on  a  high  AR  wing,  resulting  in 
a  higher  friction  drag.  Taper  ratio,  A,  would  also  show  a  greater  influence  on  the 
friction  drag  calculations,  since  on  highly  tapered  wings  much  of  the  potential  area  of 
laminar  flow  may  lie  behind  the  trailing  edge  on  the  outboard  part  of  the  wing.  Thus, 
if  Rear, t  were  to  be  increased,  the  calculated  drag  would  decrease  until  it  reached  a 
minimum  at  the  Re  which  corresponds  to  a  critical  length  equal  to  the  root  chord 
length. 

Figure  5.8  shows  the  effect  that  the  assumed  critical  Reynolds  number,  Re^n. 
hat.  on  the  calculated  friction  drag  for  the  baseline  wing.  The  friction  drag  coefficient, 
Cop  is  shown  to  decrease  as  Re„u  is  increased  to  a  value  of  about  213  million,  beyond 
which  the  upper  surfs  c.  is  fully  laminar  and  remains  oonstant.  This  is  because, 
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at  the  cruise  conditions,  an  Re^n  o i  213  million  corresponds  to  a  critical  length  of 
115  feet  (the  root  chord  length).  Using  Figure  5.8,  it  was  found  that  an  Recrit  of 
about  82  million  was  required  for  the  calculated  CD/  of  the  baseline  wing  to  be  the 
same  as  when  75%  laminar  flow  and  a  7  turbulent  wedge  were  assumed.  This  Re„.,t 
was  assumed  for  all  of  wings,  and  the  corresponding  CD/'s  were  calculated. 

Figure  5.9  shows  the  friction  drag  coefficients  as  they  were  calculated  using 
this  assumption,  along  with  those  (from  Figure  5.7)  calculated  using  the  previous 
section  s  assumptions.  It  shows  that  calculating  Co,  based  on  a  critical  Reynolds 
number  predicts  a  markedly  greater  dependence  on  both  AR  and  A.  It  may  also  be 
noted  that  the  trends  are  nearly  equal  and  opposite  to  those  predicted  assuming  a 
critical  ~  with  no  turbulent  wedge. 


5.3  Drag-Due-to-Liffc 

The  drag-due-to-lift  was  evaluated  using  the  lift  analysis  module  of  the  BDAP 
program.  This  is  the  only  drag  increment  that  is  affected  directly  by  the  camber 
of  the  wing,  thus  it  is  the  only  drag  increment  other  than  the  friction  drag  which  is 
affected  by  the  method  of  wing  design. 

Figure  5.10  illustrates  the  differences  in  the  drag  polars  of  the  baseline 
planform  cambered  using  the  two  design  methodologies.  This  figure  includes  the  wave 
drag  increment,  which  is  the  same  for  both  wings,  and  the  drag-due-to-lift  increment, 
which  of  course  vanes  with  the  camber  definition.  As  indicated  by  the  shapes  of  the 
drag  polars,  the  wing  designed  for  minimum  drag-due-to  lift  has  been  cambered  in 
such  a  way  as  to  give  it  an  advantage  at  nearly  all  positive  lift  conditions,  as  well  as 
at  the  design  lift  coefficient  of  0.08. 

The  drag-due-to-lift  increments  at  cruise  ( CL  =  .08)  for  all  of  the  wings  are 
shown  in  Figure  5.11  and  also  in  Table  5.2. 

From  Figure  5.11,  the  influence  of  taper  ratio,  A,  and  aspect  ratio,  AR,  on  the 
drag-due-to-lift  can  be  seen.  As  expected,  the  wings  with  a  higher  AR  have  a  definite 
aerodynamic  advantage.  The  HLFC  wings  are  seen  to  gain  very  little  by  increasing 
the  taper  ratio  above  about  0.15.  This  may  suggest  that  there  is  a  penalty  involved 
with  attempting  to  maintain  a  flat-roof-top  pressure  distribution  over  longer  lengths 
at  the  wing  tips. 

From  an  analysis  of  this  data,  it  was  found  that  the  differences  in  drag-due- 
to-lift  caused  by  the  design  method  are  primarily  dependent  on  taper  ratio.  At  taper 
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Table  5.2:  Drag-Due-to-Lift  at  Cruise  of  Wings  in  Parametric  Study 


Planform  ID 

Cdhlfc 

^^mirwddti 

W015 

0.00299 

0.00274 

0.9164 

W020 

0,00263 

0.00238 

0.9005 

W025 

0.00231 

0.00212 

0.9177 

W115 

0.00235 

0.00256 

0.8982 

W120 

0.00246 

0.00222 

0.9002 

W125 

0.00222 

0.00199 

0.8964 

W215 

0.00234 

0,00249 

0.8755 

W220 

0.00242 

0.00211 

0.8719 

VV225 

0.00216 

0.00190 

0.8796 

W315 

0.00287 

0.00243 

0,8467 

W32G 

0.00240 

0.00206 

0.8583 

W325 

0.00211 

0.00180 

0.8531 

ratios  of  0.0  (a  pointed  wing  tip),  the  optimized  wings  have  about  9%  less  drag  due- 
to-lift  than  the  HLFC  wings;  and,  at  taper  ratios  of  0.3,  the  optimized  wings  have 
o'-er  17%  less  drag- due- to- lift  than  their  HLFC  counterparts.  The  percent  difference 
in  drag-due-to-lift  was  found  to  have  a  nearly  linear  relationship  with  the  A  values 
between  0.0  and  0.3. 


5.4  Total  Drag 

The  total  drag,  as  explained  in  Chapter  2,  was  taken  to  be  the  sum  of  the  wave  drag, 
skin  friction  drag,  and  the  drag-due-to-lift.  Figure  5.12  shows  the  total  drag  estimates 
for  all  of  the  wings  at  the  cruise  conditions.  The  estimates  for  the  HLFC  wings  are 
shown  using  each  of  the  assumptions  of  sections  5.2.1  and  5.2.2.  It  may  be  observed 
from  these  curves  that,  in  the  analysis  of  the  total  drag,  the  differences  between  using 
a  critical  ~  and  a  critical  Re  in  the  Co,  cflculatiom*  for  the  HLFC  wings  were  greatly 
diluted. 

Of  interest  to  this  study  was  the  effect  that  aspect  ratio  and  taper  ratio  had 
on  the  drag  reductions  that  could  be  obtained  by  the  HLFC  system.  Figure  5.13 
shows  this  effect  as  a  percent  drag  reduction  as  compared  to  the  turbulent  wings 
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designed  for  minimum  drag-due*  to-lift.  Tho  results  are  shown  assuming  an  ReCri,  of 
S2  million  and  a  critical  f  of  75%  with  a  7*  turbulent  wedge.  From  this,  we  see  that 
the  assumption  of  an  Re„n  predicts  a  greater  effect  of  both  AR  and  A  on  total  drag 
reduction. 

With  the  assumption  of  a  critical  f,  it  was  found  that  the  taper  ratio  had 
little  effect  on  the  percent  drag  reduction.  In  other  words,  the  ratio  of  the  total  drag 
coefficients  was  nearly  constant  for  a  given  aspect  ratio.  However,  this  ratio  was 
found  to  be  dependent  on  the  aspect  ratio.  The  drag  reductions  range  from  slightly 
less  than  6%  to  nearly  9%  for  aspect  ratios  of  1.5  to  2.5,  respectively.  Assuming  an 
Recrll  of  32  million,  the  drag  reductions  ranged  from  slightly  less  than  5^%  up  to  8| 
-  10  %  for  aspect  ratios  of  1.5  to  2.5,  respectively.  That  is,  higher  aspect  ratio  wings 
are  likely  to  realize  greater  benefits  from  HLFC  systems  than  iow  aspect  ratio  wings. 
Thus,  they  are  more  receptive  to  such  a  design  method.  In  general,  wings  with  an 
AR  of  2.5  can  realize  up  to  twice  the  drag  reduction  through  the  use  of  an  HLFC 
system  as  a  wing  with  an  AR  of  1.5. 
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6  Concluding  Remarks 


This  study  was  undertaken  in  an  attempt  to  quantify  the  aerodynamic  benefits  and 
penalties  involved  in  designing  wings  to  accommodate  hybrid  laminar  flow  control 
(HLFC)  systems.  To  do  this,  it  was  first  necessary  to  establish  a  feasible  process 
for  the  preliminary  design  of  such  wings.  For  this  study,  an  upper  surface  pressure 
distribution  favorable  to  laminar  flow  was  investigated,  and  a  slightly  simplified  one 
was  developed  to  accommodate  the  limitations  of  the  linear  theory  used  in  the  pressure 
calculations  and  wing  design.  An  existing  design  program  was  used  to  calculate  the 
wing  lifting  load  distribution  which  would  combine  with  the  thickness  pressures  to 
produce  this  desired  upper  surface  pressure  profile  while  maintaining  a  spanwise  lifting 
distribution  designed  to  minimize  drag-due-tc  lift.  Linear  theory  was  then  used  to 
define  the  camber  definition  required  for  the  lifting  load  distribution. 

This  design  methodology  and  a  method  used  to  camber  wings  for  minimum 
drag-due-to-lift  were  applied  to  several  planform  shapes  in  an  attempt  to  evaluate 
the  effects  of  aspect  ratio  and  taper  ratio  on  the  potential  performance  benefits  of 
HLFC  wings.  The  skin  friction  drag  was  shown  to  be  almost  independent  of  AR  and 
A  for  the  HLFC  wings  at  about  30%  lower  values  than  those  for  the  turbulent  wings, 
if  it  were  assumed  that  75%  of  the  chord  length  on  the  upper  surface  were  laminar 
except  for  a  7°  wedge  at  the  root.  The  friction  drag  was  also  calculated  for  all  of  the 
wings  assuming  a  critical  Reynolds  number  of  82  million.  With  this  assumption,  Cp, 
was  predicted  to  decrease  with  increases  in  either  AR  or  A.  The  differences  in  Cp, 
between  the  HLFC  wings  and  their  turbulent  counterparts  were  also  predicted  to  be 
much  greater  for  the  high  AR  and  high  A  wings. 

It  was  also  calculated  that  the  wings  designed  for  HLFC  had  9-17j%  more 
drag-due-to-lift  than  the  corresponding  wings  designed  for  minimum  drag-due-to- 
lift.  This  percentage  was  found  to  linearly  increase  with  A  and  to  be  relatively 
independent  of  AR.  However,  when  this  drag  increment  was  combined  with  the 
wave  drag  and  friction  drag  increments,  the  tota)  drag  benefit  of  the  HLFC  wings 
was  mostly  dependent  upon  the  wing  AR.  This  is  because  the  unfavorable  effects  on 
diag-due-to-lift  increase  with  A  at  about  the  same  rate  as  the  beneficial  effects  on  the 
skin  friction  drag  (assuming  75%  laminar  flow)  for  the  HLFC  wings. 
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In  conclusion,  it  is  predicted  that  the  total  drag  reductions  accomplished  by 
designing  wings  for  HLFC.  as  opposed  to  designing  them  for  minimum  drag-due-to- 
lift,  vary  from  5|%  to  nearly  9%  based  on  a  critical  -  of  75%  and  from  5%  to  nearly 
10%  based  on  /?etr(,=82  million  depending  on  the  wing  A  and  AR.  Experimental 
investigation  is  required  to  resolve  which  of  these  assumptions  pertaining  to  ihe 
location  of  transition  to  turbulence  is  more  correct.  It  was  also  found  that  this 
relationship  was  primarily  dependent  on  AR  and,  in  general,  wings  with  higher  aspect 
ratios  are  likely  to  realize  greater  benefits  from  HLFC  systems  than  low  aspect  ratio 
wings. 

The  next  logical  step  in  this  study  would  be  to  implement  higher  order  methods 
in  the  inverse  design  process.  But  until  such  methods  are  made  available,  progress 
must  be  made  with  the  tools  at  hand.  The  present  research  is  not  intended  as  a 
substitute  for  higher  order  methods;  it  is  hoped  instead  that  it  will  augment  interest 
in  the  advancement  of  such  methods  so  that  they  may  be  used  in  future  studies  and 
to  verify  the  present  results. 


to 
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Drag 


Figure  2.1  -  Superposition  Method  of  Drag  Analysis 
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(a)  series  of  Mach  plane  cuts  used  in  the  calculation  of  a  single  equivalent  body 


(b)  equivalent-body  area  distribution  for  a  sample  0  ~  0\ 


Figure  2.2  -  Far-Field  Wave  Drag  Calculations 
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Figure  2.3  -  Mach  Box  Computational  Grid  for  Supersonic  Aerodynamics 
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*Kf  is  a  turctioii  of : 

-Airfoil  geometry 
•  Normal  Mach  number 
-Reynolds  number 


•  Attainable  thrust  factor,  Ky  =  T! 
Kf  =  1.0  —  Attached  Sow 
Kr  =  0.0  —  Vortex  flow 


en  =  — - — sm(cos_1/f,7'  —  <5) 
coi  Afc 

Ca  —  -C(  cos(coj“l  Kf  -  S) 

Figure  2.4  -  Theory  of  Attainable  Vortex  Thrust 


Figure  2.5  -  Mixed  Boundary  Layer  Schematic 
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Figure  2.8  -  Continued 
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Figure  3.1  -  NACA  64-004  Airfoil  Thickness  Profile 
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Figure  3.3  -  Pressure  Distribution  on  Baseline  Wing  Designed  for  Minimum  Drag 
Due- to- Lift 


4.1  -  Lifting  Pressure  Distribution  on  Flat-Roof- Top  Camber  Surface 

(combines  with  thickness  pressures  to  give  a  flat-roof-top  Cp  distribution) 


Figc  ■*  4.2  -  Lifting  Pressure  Distribution  on  Span-Optimized  Camber  Surface 
(optimized  in  the  presence  of  the  loading  of  figure  4.1) 


(a)  Flat-Roof-Top  Camber  Surface 


(b)  Spanwise-Optimized  Camber  Surface 
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Figure  4.4  -  Pressure  Distribution  on  Uncambered  Baseline  Wing  at  Design  Lift 
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Figure  4.11  -  Pressure  Distribution  on  Baseline  HLFC  Wing  with  Root  Section 
Rotated 
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Figure  5.4  -  Effect  of  Laminar  Flow  on  FYiction  Drag  of  Baseline  Wing 
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Hgurc  o.5  -  Effect  of  AR  and  taper  ratio  on  Skin  Friction  Drag  of  Fully  Turbulent 
Wing 


Machs=2.0 


Figure  5.6  -  Effect  of  AR  and  taper  ratio  on  Skin  FVictione  Drag  of  Fully  Laminar 
Wing 
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Figure  5.7  - 


MachsZO 


Effect  of  Turbulent  Wedge  at  Root  on  Friction  Drag 
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figure  5.8  -  Effect  of  Critical  Reynolds  Number  on  friction  Drag  of  Baseline  Wing 
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Figure  5.9  -  Effect  of  Critical  Reynolds  Number  Assumption  on  Calculated  Skin 
Friction  Drag 
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Figure  5.10  -  Near-Field  Wave  Drag  Plus  Drag-Due- to-Lift  on  Baseline  Planform 
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Figure  5-11 


Cruise  Drag-Due-to-Uft  of  Parametric  Wings 
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Figure  5.12  -  Total  Drag  of  Parametric  Wings 
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Figure  5.13  -  Total  Drag  Benefit  of  HLFC  Wings 
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